Nonlinear photonic lattices in anisotropic nonlocal self-focusing media 
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We analyze theoretically and generate experimentally two-dimensional nonlinear periodic lattices 
in a photorefractive medium. We demonstrate that the light-induced periodically modulated non- 
linear refractive index is highly anisotropic and nonlocal, and it depends on the lattice orientation 
relative to the crystal axis. We discuss stability of such induced photonic structures and their 
guiding properties. 



The study of nonlinear effects in periodic photonic 
structures recently attracted strong interest because of 
many novel possibilities to control light propagation, 
steering and trapping. Periodic modulation of the refrac- 
tive index modifies the linear spectrum and wave diffrac- 
tion and consequently strongly affects the nonlinear prop- 
agation and localization of light |l| . 

Photonic lattices can be optically induced by linear 
diffraction-free light patterns created by interfering sev- 
eral plane waves 0. However, the induced change of 
the refractive index depends on the light intensity and, 
in the nonlinear regime, it is accompanied by the self- 
action effect [2j. The nonlinear diffraction- free light pat- 
terns in the form of stable self-trapped periodic waves 
can propagate without change in their profile, becoming 
the eigenmodes of the self-induced periodic potentials. 
This behavior is generic, since nonlinear periodic waves 
can exist in many types of nonlinear systems, and they 
provide a simple realization of nonlinear photonic crys- 
tals. Such structures are "flexible" because the lattice 
is modified and shaped by the nonlinear medium; these 
flexible lattices extend the concept of optically-induced 
gratings beyond the limits of weak material nonlinear- 
ity. Moreover, the nonlinear lattices offer many novel 
possibilities for the study of nonlinear effects in periodic 
systems because they can interact with localized signal 
beams via the cross-phase modulation and can form a 
composite bound states 0, || . 

Nonlinear photonic lattices created by two-dimensional 
arrays of in-phase solitons have recently been demon- 
strated experimentally in parametric processes H, and 
in photorefractive crystals with both coherent |7| and 
partially incoherent 0-0-0 light. For the case of two- 
dimensional arrays of in-phase solitons created by the 
amplitude modulation, every pixel of the lattice induces 
a waveguide which can be manipulated by an external 
steering beam PUSH- However, the spatial periodicity 
of these lattices is limited by attractive soliton interaction 



that may leads to their strong instability. In contrast, 
the recently suggested two-dimensional lattices of out-of- 
phase solitons are known to be robust in isotropic sat- 
urable model ^3 ■ The phase profile of such self-trapped 
waves resembles chessboard with the lines of 7r-phase 
jumps between neighboring white and black sites. 

In this Letter, we study two-dimensional nonlinear lat- 
tices with a chessboard phase structure in anisotropic 
nonlocal self-focusing media and generate such lattices 
experimentally in a photorefractive crystal. We demon- 
strate that the light-induced periodically modulated non- 
linear refractive index is highly anisotropic and nonlocal, 
and it depends on the lattice orientation relative to the 
crystal axis. We discuss stability of such induced pho- 
tonic structures and their guiding properties. An obvious 
advantage of using this novel type of nonlinear periodic 
lattices when compared with in-phase lattices or incoher- 
ent soliton arrays 0,0 is that such lattices can be made 
robust even with smaller lattice spacing. This creates 
an opportunity for experimental studies of mutual cou- 
pling of periodic and localized waves recently discussed 
for one-dimensional geometry, as well as the generation 
of two-dimensional composite lattice solitons 

Spatially periodic nonlinear modes appear naturally 
due to self-focusing effect and modulational instabil- 
ity 0. When self- focusing compensates the diffraction 
of optical beams, it may support both isolated spatial 
solitons and periodic soliton trains, or stationary peri- 
odic nonlinear waves. The latter include well studied 
cnoidal waves, described by the Jacobi elliptic functions 
en and dn, solutions to the nonlinear Schrodingcr (NLS) 
equation, 



id z E + V 2 E + n(I)E = 0, 



(1) 



where / = \E\ 2 and, for one-dimensional case, V 2 = 
d 2 /dx 2 — d 2 . Similar stable periodic waves exist in 
different nonlinear models, including quadratic, Kerr- 
type saturable nonlinearities etc. The family of two- 



dimensional nonlinear periodic waves | 10| can be also ex- 
tended to the case of the rectangular geometry with two 
different transverse periods, because such anisotropic de- 
formations of the square lattice do not enhance its modu- 
lational instability. Stabilization of phase engineered soli- 
ton arrays was reported recently for anisotropic model @ . 

In this paper, we consider photorefractive crystal as 
an example of anisotropic and nonlocal nonlinear media. 
In this case, the nonlinear contribution to the refractive 
index in the model is given by 

n(i) = ra^, (2) 

where the electrostatic potential ip of the optically- 
induced space-charge field satisfies a separate equation: 

Vj_9? + Vx</>Vxln(l + J) = d x hx(l + I). (3) 

Here is the two-dimensional gradient operator, = 
d 2 . + d 2 ., and the intensity / is measured in units of the 
background illumination (dark) intensity, necessary for 
the formation of spatial solitons in such a medium. The 
physical variables i, y, and z correspond to their dimen- 
sionless counterparts as (x, y) — xq(x, y) and z = 2kxqZ, 
here x$ is the transverse scale factor and k — 27rno/A 
is the carrier wave vector with the linear refractive in- 
dex uq. Parameter T = x\n?"n^r c fiE is defined through 
the effective electro-optic coefficient r c g and externally 
applied bias electrostatic field £. 

Stationary solutions to the system CQ)-© are sought 
in the standard form, E(x, y,z) — U (x, y) exp(ifcz) where 
the real envelope U satisfies the equation 

-kU + V\U + Td x ipU = 0. (4) 

We look for periodic solutions, U(x, y) = U(x + 2w, y + 
2tt), and solve Eqs. ©, using the relaxation tech- 
nique 1 1 1 1 ] with initial ansatz in the form of a linear pe- 
riodic mode, Un n (x,y) — v4sin:r siny. We find that at 
least two distinct families bifurcate from the linear wave 
Unn(x,y), depending on its orientation: a square pat- 
tern parallel to the c-axis, and a diamond pattern ori- 
ented diagonally. Figures 1(a) and 1(b) show the field 
and refractive index distributions in the low (fc = —1.9, 
A ~ 0.9) and relatively high (k = —1.5, A ~ 3.6) satura- 
tion regimes for both families. In a general case r ^ 1, 
these two families occupy a band k 6 [— 2, T — 2] with 
the amplitudes A(k) and power densities P(k) vanish- 
ing in the linear limit k — > —2, see Fig. 1(c). Here 
the power density is defined as the power of a unit cell, 
P = Ajf* U 2 dxdy. The main difference between two 
solutions, clearly seen in Figs. 1(a) and 1(b), comes from 
the refractive index: the regions with the effective fo- 
cusing lenses are well separated for the diamonds, and 
fuse to the effectively one-dimensional stripes for square 
pattern, in the limit of strong saturation. In Fig. 1(c), 
we plot maximal and minimal values (Extrema) of the 
refractive index, Extr^^y). 
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FIG. 1: (a) Square and (b) diamond self-trapped stationary 
periodic patterns in the model flT)-J3t at T = 1. Families of 
nonlinear waves are summarized in (c) and (d) ; solid (dashed) 
lines correspond to the square (diamond) patterns. Both lines 
for A, P, and max(d m <p) coincide in (c). Horizontal dotted 
lines in (d) correspond to the linear limit F — > 0. 

In Fig. 1(d) we show the FWHM of a single intensity 
spot in two orthogonal directions, d x ,y, characterizing the 
degree of spatial localization of cnoidal wave ^(|- The 
ellipticity of every lattice site, e = d y /d x > 1, depends 
on the propagation constant k, similar to the ellipticity 
of a single photorefractive soliton see Fig. 1(d). 

In order to test the lattice stability, we propagate nu- 
merically two types of initially perturbed periodic solu- 
tions and observe robust propagation for the distances ex- 
ceeding the experimental crystal length. Figure 2 demon- 
strates an example of a stable propagation for the param- 
eters close to our experimental situation. 

To demonstrate experimentally both existence and sta- 
bility of these nonlinear periodic lattices in anisotropic 
and nonlocal media, we use a setup similar to that 
employed earlier [f|. Linearly polarized beam from a 
frequency-doubled Nd:YAG laser at 532 nm is sent to 
a liquid crystal programmable spatial light modulator in 
order to create a periodic light pattern with a variable 
period and orientation. The output of the modulator 
is then imaged by a high numerical aperture telescope 
(demagnification ~ 10) on the front face of a photore- 
fractive SBN:60 crystal. The incident light is linearly 
polarized parallel to the c-axis, thus experiences strong 
photorefractive nonlincarity. The imposed pure phase 
modulation transforms into an amplitude modulation of 
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FIG. 2: Numerical results for the propagation of (a) square 
and (b) diamond self-trapped patterns for T = 11.8 in the 
low saturation regime A ~ 1, k — —0.5. Intensities of the 
lattice and probe wave it guides are shown after propagation 
z ~ 23 mm. On the input, the lattice is perturbed with 20% 
of random noise, and the probe is a broad Gaussian beam. 
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FIG. 3: Experimental results for the (a) square and (b) dia- 
mond lattices. Left to right: output intensities of the lattice 
after linear and nonlinear propagation, and output of a guided 
plane wave, respectively. 

the beam at the front face of the crystal, where noise 
is reduced by proper spatial filtering. The crystal is ex- 
ternally biased and uniformly illuminated with a white- 
light to control the dark irradiance. The generated peri- 
odic pattern experiences robust linear propagation in the 
crystal (at zero bias field) with a negligible change in the 



periodicity. The linear output of the generated periodic 
wave with periodicity of 31/im at the back face of the 
crystal is shown in Fig. 3 (left column) for two different 
orientations with respect to the crystal axis. 

Applying an external DC electric field (£ = 1000 
V/cm) across the crystal creates appropriate conditions 
for the formation of spatial solitons and, correspondingly, 
influences the propagation of the periodic waves. Since 
with increasing of the nonlinearity (i.e., biasing voltage) 
the output of the lattice does not change significantly 
[Fig. 3 (middle column)], the only way to test that the 
lattice propagates nonlincarly is to probe the periodic 
potential it induces. To do so, we send a broad plane 
wave trough the crystal and observe its modulation at 
the output. In practice, this is realized by switching-off 
the voltage on the modulator, thus removing the mod- 
ulation imposed on the light pattern and generating a 
broad plane-wave at the input of the crystal. This ap- 
proach assures that the plane wave is propagating exactly 
along the induced waveguides. Due to a slow response of 
the photorefractive crystal, we can quickly monitor the 
output of the plane wave without modifying the induced 
refractive index change. Output intensity distributions 
for two orientations of the lattice pattern are shown in 
Fig. 3 (right column) . The plane- wave gets split into sev- 
eral channels, and we observe guiding of the probe beam 
at the maxima of the refractive index. The experimental 
pictures shown in Fig. 3 are in a good agreement with 
the corresponding numerical simulations of the nonlin- 
ear anisotropic nonlocal model presented in Fig. 2, and 
demonstrate a qualitative difference of the guided pat- 
terns for two different orientations of the lattice. 

In conclusion, we have studied theoretically and gener- 
ated experimentally two-dimensional nonlinear photonic 
lattices in anisotropic photorefractive medium. We have 
found two distinct classes of self-trapped robust spatially- 
periodic waves with out-of-phase neighboring sites, the 
square pattern oriented parallel to the crystal axes, and 
the diamond pattern oriented diagonally in the transverse 
plane. We have demonstrated that the highly anisotropic 
refractive index distribution induced by the lattice differs 
significantly from its isotropic counterpart and depends 
strongly on the lattice orientation. 
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